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Abstract 

Laser interferometer gravitational wave detectors can be operated at their 
free spectral range frequency. We show that in this case and when the inter- 
ferometer is weh understood one could detect a stochastic background using a 
single detector. 

1. Introduction 

We define the stochastic gravitational wave background (sgwb) by a dimensionless 
strain h(t) which is a real random variable, stationary in time. Namely, it is Gaussian 
distributed with zero mean 

< h{t) > = (1) 

and finite "power" 

< h^{t) > = P (2) 

Eqs.(l,2) do not suffice to completely describe the random process. To do so we must 
also specify how the power P is distributed in frequency^, what is referred to as the 
power spectral density (PSD). 

We therefore assume that h{t) can be represented, at a given point in space, as a 
Fourier integral 

/oo 
h{f)e-^-^f'df (3) 
-oo 

Our assumptions about the properties of h{t) imply that 

< kf) > = (4) 

and that [1] 

< h\f)h{f') > = ^-6{f - f')S{f) < / < oo (5) 



^For instance whether the random process is like "white noise", "1/f noise", resonant at a given 
frequency, etc. 
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The function S{f) is the (one-sided) power spectral density per unit time. The proof 
of Eq.(5) and and the subtleties in the definition of h{f) are discussed in the Appendix. 
Throughout this note we have adopted the nomenclature and conventions of reference 
[2] as much as possible. 

The 5-function can be removed from Eq.(5) by integrating over /' in the narrow 
interval / - A//2 < /' < / + A//2, to obtain 

/-/+A//2 ^ ^ 

-^a) =A/-.o 2 / ^ <h*{f)h{f)>df^2<\h{f)\'>Af (6) 

J/-A//2 

The integration over the small interval A/ expresses a kind of smoothing of h{f) [3,4]. 
If h{f) is obtained from a stretch of data of length T, then A/ = 1/T and Eq.(6) 
takes the form 

Sif) = ^<\Hf)\'> 0</<oo (7) 

It follows from Parceval's theorem that 



/ S{f)df = -/ < \h{f)\' > df ^ - \h{t)\'dt^ 

Jo 1 J-oo 1 J-T/2 



(8) 



= <\h{t)\'>^P 



Eqs.(7,8) are our central result. Experimentally we wish to determine S{f). We will 
introduce the notation 



kf) = ^^sU) (9) 

and refer to h{f) as the sgwb amplitude (density) per square root of frequency, i.e. 
h{f) represents strain /\/Hz. 



2. Single noiseless detector 

In an experiment we always deal with finite time intervals and with time series 
(rather than functions of time). This is well suited to the form introduced in Eqs.(6,7). 
We consider a linear detector that in response to g.w. strain h{t) outputs a signal 
a{t). To include the response of the detector we must work in the frequency domain 

A{f) = H{f)h{f) (10) 
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where H{f) is a deterministic transfer function^ peaked at some value /o and decreas- 
ing below and above /q. The detector output is given by 



/oo 
Aif)e-'^'f'df (11) 
-oo 

The output a{t) is sampled at a rate 1/A = 2/^, where fc is referred to as the 
Nyquist critical frequency. We examine a record consisting of N data points; the 
length of the record is r = A^A and the resulting time series is 

an = a{tn) with t„ = nA n = 0, 1, . . . , (A^ - 1) (12) 

We can form the Discrete Fourier Transform (DFT) of this time series [2] . 

N-l 

A = E Qne^-'^"/^ 0,1,..., (N-l) (13) 

n=0 

This is a mapping of the complex numbers a„ (of course in our application the 
are real) onto N complex numbers A/-. The mapping is invertible 



1 Af-l 

a„.^4.Yl ^.e-^-^"/^ (14) 



and obeys Parceval's theorem 



N-l -1 N-l 

n=0 k=0 

The elements Ak are related to A{fk), the Fourier transform of a{t), through 

A{fk) ^ AAk (16) 

We can obtain an estimate of the PSD by forming the periodogram of a{t). The 
frequency is defined only for zero and for N/2 positive values 

/fc = ^=2/,A k = 0,l,...,N/2 (17) 

The bandwidth B is the frequency spacing between k and A; + 1, or 

B = 1/NA = 1/t (18) 

^We take H{f) to be dimensionless, that is the transfer function from g.w. strain to strain recorded 
by a cahbrated detector. 
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The terms of the periodogram are 



P(0) 



7V2 



A; = l,2,...,(7V/2-l) (19) 



^ \A P 



For a{t) real = l^-fcP and therefore 



(20) 



We now use the approximate form of Eq.(16) in Eq.(7) to obtain for the PSD per 
unit time 



2 2A2 
SAifk) = -<\A{h)\'>^ <\A,\'> 

T T 



- r^.<\M'> = ip{h) 



(21) 



^2 I «l 5 

Returning to Eq.(lO) we can write for the discrete case 

AAk ^ A{fk) = H{h)h{fk) (22) 

and therefore 

SAifk) ^ -\H{h)? <\h{fk)?> (23) 

T 

If h{f) is approximately constant in the region where H{f) peaks then the detected 
spectral density (per unit time) SA{fk) will follow the spectrum o/|iJ(/fc)|^. 



3. Single noisy detector 

Suppose now that the detector introduces noise so that the output signal is 

c{t) = q{t) + a{t) (24) 

Here q{t) is a stochastic variable with zero mean and standard deviation a. After 
sampling c{t) and forming the DFT we obtain as in Eq.(13). 
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Ck — Qk-\- -^k 



(25) 



Since q{t) and a{t) are uncorrelated, the PSD/time, upon averaging, will contain only 
two terms 



ScUk) - 



< \Qk? > + < 



> 



(26) 



If the noise spectrum is flat in frequency, one can extract < |^fcp > by fit- 
ting Sc{fk) with the expected spectral shape, which according to Eq.(23) is that of 
|H(/fc)p. In this way one can determine < \h{fk)\^ >, in the region fk — /o, with a 
single detector, provided the spectrum of < IQfcP > is known. 

In fact, one can work directly with the magnitude of Ak which is related to /i(/), 
defined by Eq.(9). To extract a meaningful signal, \Ak\ must be larger or equal 
than the fluctuations of Qk (rather than Qk itself). For a single measurement the 
fiuctuations of \Qk\ are equal to < IQfeP > since q{t) was assumed to have zero 
mean. After averaging Na different DFT spectra 



1 



■\Qk 



(27) 



Further, in fitting the data to extract \Ak\'^ we have at our disposal the n frequency 
bins in the bandwidth A/, typically a few times the FWHM of H{f), 



n = Af/B = tA/ (28) 
The significance of the fit can be expressed as a signal to noise ratio 



N 



\Na - 



J< \Ak\'^ > ^ I — I 

J< \Qk\' > 



< \Ai 



> 



<\Qk\'> 



y/TAf\Hk 



I'M. 



(29) 



In Eq.(29) we have set T = rNa for the total measurement time and \Hk\ is the 
magnitude of the transfer function averaged over the frequency band A/. We also 
used Eqs.(25,26) and the definition of Eq.(9) to relate < l^lfel > and < \Qk\ > to 
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the amplitude (densities) per unit frequency for the sgwb signal and the noise in the 
detector. 

Eq.(29) shows that when using a single interferometer, the S/N ratio for hg/hN 
improves as ^/TAf, in contrast to the correlation method [see Eq.(37)] where the 
S/N for the same ratio improves only as (TA/)-*^/^. This is particularly true if no 
signal is detected and the measurement yields only an upper limit on hs- Computer 
simulations [5] of a sgwb signal using for noise the actual time series from LIGO test 
data confirm this conclusion. 



4. Correlated detectors 

The "standard" approach to the measurement of a stochastic signal is to correlate 
the output of two detectors. For the sgwb this method has been discussed extensively 
in [6,7] and has been applied to data from the LIGO interferometers in [8]. Here we 
only briefly outline the analysis procedure and present the conclusions. 

Assuming that the two detectors 1,2 are co- located and co-ahgned, the stochastic 
signal, a{t) will be the same in both detectors whereas the noise signals qi{t) and 
q2it) are treated as uncorrelated. Thus we write 



ci{t) = qi{t) + a{t) 

(30) 

C2{t) = q2{t) + a{t) 

We form the cross-correlation (at zero lag) of the two signals and integrate over a 
time interval r, to obtain a statistic [1] 

r ci{t)c2{t)dt (31) 
Jo 

The mean value /i, of this statistic obtained from repeated cross-correlation spectra 
(samples) is a measure of the stochastic signal because 

<C> = r < a^it) >dt^T < a^it) > (32) 
Jo 

We can express in terms of frequency domain variables by using the definition of 
Eq.(8) 
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fOO 

li^T SA{f)df 

Jo 



(33) 



where 5'^(/) is the spectral density of the stochastic signaL 
The fluctuations in n are given by the variance of C 




(34) 



Jo 



/ < ql{t) >< ql{t) >dT^T 



/ S^,{f)S^,{f)df 
Jo 



(35) 



where Smif) are the spectral densities of the noise. If we use Ns samples to determine 
the mean, then < /i > will be Gaussian distributed around its true value with standard 
deviation ac / V^s- Thus the signal to noise ratio is 



If we further assume that 5'a(/) is reasonably constant in the region of integration, 
and so are <S'jVi(/) — S^^lf) we can use the notation of Eq.(9) to simplify Eq.(36) 



where wc set T = Nst for the total measurement time. Eq.(37) is to be compared 
with Eq.(29) obtained for measurements with a single detector. 
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APPENDIX A 



Spectral Density of a Random Variable 

We expressed a stochastic signal h{t) by Eq.(3) 

/oo _ 
h{f)e-'^'f'df (Al) 
-oo 

It is not possible to invert Eq.(Al) because the integral 

/oo 
h{t)e^^'f'dt (A2) 
-oo 

diverges for a stationary signal. 

Instead, following [4] we can introduce the truncated function 

hrit) = h{t) -T/2<t< T/2 

— elsewhere (A3) 



Thus 



/oo 
hT{t)e^^'^'dt (A4) 
-oo 



and we form 



h*T{f)hT{f) = / / h*T{t)hT{t')e-^'''^'e^'''^' dtdt' (A5) 

J —oo J —oo 

We seit' — t-\-T so that 

h*T{f)hT{f') = / / h*T{t)hT{t + r)e^^'^f~f^'e^^'f^dtdr (A6) 

J —oo J —oo 

We now take the ensemble average of both sides, interchange the order of averaging 
and integration on the right-hand side, and define the autocorrelation function [9, 3] 

R^{t) = < h*T{t)hT{t + r) > (A7) 
The function Rt{t) vanishes for |t| > T, and Eq. (A6) reads 

~ ~ poo 

< hUmAf) >^S{f- f) / RT{T)e''^'f^dT (A8) 

^— oo 
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We therefore define the power spectral density (per unit time) as the Fourier transform 
of the autocorrelation function 

/oo 
RT{r)e'^'''^^dT < / < oo (A9) 

-oo 

Thus Eq.(A8) is recast as 

< hUDhAf) > = l^if - f')ST{f) (AlO) 

By letting T — > oo we obtain Eq.(5) of the main text. 

To prove Eq.(7) of the main text we proceed as follows. By definition, [Eq.(A7)], 

Rt{0) = < \hT{t)\' >^7f; l \hTit)\'dt = - / \hT{t)\'dt = 

1 J-T/2 1 J-oo 

-Tl_J^T{f)?df (All) 
Integrating Eq.(A9) over frequency we obtain 



/■oo 1 fCO nOO 

/ Srim = - / 2RT{T)^^'^''dfdT = / RT{T)5{T)dT = i?T(0) = 
JO ^ J— oo J— oo 

where we used Eq.(All) for -Rr(O). Equating the integrands we establish 

Srif) = ^\hT{f)\' (A13) 

which is analogous to Eq.(7) of the text but for the Fourier transforms hxif) of the 
truncated funcfAons hT{t). It woTild seem that by letting T — oo we could recover the 
relation shown in Eq.(7). This is not the case [10] because Srif) does not converge 
to a single value as T — > oo. It can however be "smoothed" by taking the ensemble 
average of |/iT(/)P before letting T — > oo, so that 

S{f)='r-.ool<\hT{f)?> (A14) 
which is the mathematically precise statement of Eq.(7). 



9 



It is of some interest to express Eq.(AlO) in the time domain. Using the inverse 
of Eq.(A4) we write 

/ h\f)e'^'H{r)e--'^'f'Ufdf (A15) 
Ensemble average, interchange the order of integration and averaging, and use Eq.(AlO) 

1 /' oc f oo 

< h*^{t)hT{t') > = - / 6{f- f)ST{f)e'''^^f'-^'''Ufdf 

Z J— oo J —oo 

1 pOO 

= ^ / e'^'f^'-'">ST{f)df (A16) 

Z J —oo 

When t' — t, Eq.(Al6) reduces to Eq.(8) of the main text in the hmit T — > oo. 

The physical interpretation of the autocorrelation function, Eq. ( A7) . implies that 
if the signal h{t) "repeats" itself, on the average, with period Tq, then it must contain 
a frequency component / = I/tq. For example if h{t) — coscuot — cos(27r/oi), then 

Rt{t) = cos(27r/oT) and Srif) = 6{f - fo) (A17) 

If Rt{t) drops off for r > tq, then the highest frequency contained m.h{t) is / < I/tq. 
For white noise we can write Rt{t) — Tc5(t) where Tc is the correlation time, defined 
for instance by our sampling interval. Then the spectral density is flat up to a 
frequency /c = l/n. 
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